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The 7/ signal from partially quenched Wilson fermions 
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We present new results from our ongoing study of flavor singlet pseudoscalar mesons in QCD. Our approach 
is based on (a) performing truncated eigenmode expansions for the hairpin diagram and (b) incorporating the 
ground state contribution for the connected meson propagator. First, we explain how the computations can be 
substantially improved by even-odd preconditioning. We extend previous results on early mass plateauing in the 
r] channel of two-flavor full QCD with degenerate sea and valence quarks to the partially quenched situation. We 
find that early mass plateau formation persists in the partially quenched situation. 
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1. Introduction 

Hadron spectroscopy has been the testbed for 
monitoring the progress of lattice QCD calcula- 
tions for over more than two decades. The rj me- 
son has been rather evasive, however pl^P^,||. 
Being a flavor singlet object, its correlation func- 
tion is a combination of a connected (octet) con- 
tribution and the infamous hairpin diagram. The 
latter is difficult to evaluate, for it is the correla- 
tor of two trace expressions over the inverse her- 
mitean Dirac operator, Q^ 1 . 

So far we have been computing D by means of 
a truncated eigenmode representation of Q^ 1 , 
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with Q\i/)i) = and N = dim(Q). 'Trun- 

cated' means that we restrict the sum Eq. (|l|) to 
i -C N, i.e. to some low lying eigenmodes 0. An- 
other important ingredient of our approach is a 
ground state projection of the connected corre- 
lator (one loop ground state analysis = OLGA), 
prior to the combination with the data for the 
hairpin-operator. 

For the 16 3 x 32 SESAM lattice at f3 = 5.6 we 
showed for 2 different quark mass values that the 



300 lowest lying eigenmodes of Q suffice to ap- 
proximate D well enough. With OLGA, we found 
long rj' mass plateaus with onset at the very first 
time slice. The results appear to be fully con- 
sistent with former stochastic estimator analyses, 
which provides evidence that truncation effects 
from Eq. ([!]) in the SESAM operating conditions 
Q can be safely discarded. 

We observed in Ref. [[!] that stochastic esti- 
mator techniques and TEA require comparable 
computational effort. Actually, we can reduce the 
costs of TEA by means of preconditioning. 

2. Even-odd preconditioning of Q 

Since the Wilson Dirac matrix M connects 
nearest neighbor sites only, it can be readily even- 
odd preconditioned, i.e. cast into the form 



M = 
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where e(o) stands for even(odd). Note that 75(1 — 
n 2 D oe D eo ) is a hermitean matrix with orthogonal 
eigenvectors 4>i and real eigenvalues <Tj. In terms 
of these eigenmodes, Q^ 1 = (j^M)^ 1 reduces to 
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For brevity we omitted a traceless term 75 l c 
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Eq. (Q) now becomes 
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TtQ- 1 =TrJ2- 
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The quark loop signals as determined from the 
two different sets of vectors, i.e. from Eq. (|l|) and 
Eq. ([!]), are practically indistinguishable. Note 
that M offers two crucial benefits: 1. it requires 
only half the memory and 2. the eigenmode algo- 
rithm applied on 75(1 — K 2 D oe D eo ) is improved 
by a factor 4 in convergence rate. In fact the or- 
der of the Chebyshev (accelerating) polynomials 
can be cut down from 80 to 20, see ||. 

3. Unquenched 2-fiavor rf mass analysis 

Let us present next an update of our recent 
unquenched 2-flavor ij' mass analysis 0j . We have 
by now completed the quark loop computations 
at 4 different sea quark mass values, each gauge 
field ensemble containing 0(200) gauge fields. 

The local 7/ masses as extracted from logarith- 
mic derivatives of the ?/-correlator are exhibited 
in Fig. ([I]). We find that OLGA produces long 
and stable mass plateaus. As expected the quality 
of the plateaus improves when going more chiral 
in sea quark masses. 

In Fig. (^) we plot the chiral extrapolations for 
the singlet masses. Our data clearly favor the 
quadratic fit, which yields \ 2 — 0.14, compared 
to x 2 = 5.5 for the linear extrapolation. 

We retrieve a definitely non vanishing mass 
gap between singlet and nonsinglet masses after 
quadratic chiral extrapolation. The value of the 
computed mass gap is far from the experimental 
one (even if we scale the Ty'-mass according to the 
Witten-Veneziano expectation down to a pseu- 
doexperimental value as indicated in Fig. (|^)). 
Physically this could be due either to limitations 
of unimproved Wilson fermions or the omission of 
strange quarks. Here we address the latter. 

4. Partially quenched rj — r)' mass analysis 

Given the SESAM gauge fields with just two 
dynamical flavors, let us turn now to a partially 
quenched scenario (PQS) we introduce 3 

quenched quarks, that will be referred to as u, 
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Figure 1. r\' mass plateaus from unquenched 
QCD with two degenerate flavors, plotted in lat- 
tice units, for 4 values of n sea . The lines corre- 
spond to the octet masses. 
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Linear chiral extrapolation, 2 flavours 
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Figure 2. Chiral extrapolations of the pseu- 
doscalar octet and singlet masses in unquenched 
QCD. 

d and s, with m u — m^. For u and d we will 
choose K va i , independently from n sea . In this set- 
ting we can extrapolate eventually to light u and 
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Figure 3. Local rj masses in the u-d sector, 
Eq. (Q), from PQS, for 3 values of n va i- The lines 
correspond to the octet masses. 
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Figure 4. Same as Fig. (||) but for the s sector. 



the sea quark mass of the QCD vacuum configu- 
rations. 

For any value of time separation t, the mass 
gap, mo(t), can be extracted by fitting the mo- 
mentum zero Fourier transform of Eq. (^|) to 
the lattice data, once the local nonsinglet pseu- 
doscalar meson masses rrii(t) have been deter- 
mined from a standard lattice analysis. The rf 
masses in the u — d, s sectors follow subsequently 
from Eq. (Q), Eq. (g), respectively. 

It is nontrivial to trace the rf signals under par- 
tial quenching, since we see statistically well de- 
fined linear and constant terms in the propaga- 
tors, reflecting the double pole structure. Their 
effect increase as we move m va i away from m sea . 
To illustrate the actual rf signal from PQS, we 
plotted the mass plateaus of rf = (uu + dd) 
for K S ea = 0.1575 and three different n va i val- 
ues in Fig. (|^). Increasing m va i, the plateau 
signal remains stable. Similarly, Fig. (^) shows 
the plateauing in the s sector as computed from 
Eq. (|). 

The quality of our plateaus gives us confidence 
that we will be able to successfully complete the 
partially quenched calculations and move on to 
the study of improved fermions. 
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d valence quarks at any given value of n sea . 

In PQS, the momentum space Green's func- 
tions, corresponding to the states uu, dd and ss, 
are given by 



Gij — 



2ml 
p 2 + rnj- 
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(6) 
(7) 

*o - ,/ v - TO? , T)' = {ss) . (8) 

In these formulae, the indices i,j refer to the va- 
lence quark masses, while the index / points to 



p 2 + TO? [p 2 + TO?) {p 2 + TO?) F(p 2 

with 

F(p 2 ) = 1 + 
and 



2ml — m ff ~ m h v' = ( uu + dd) 
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